Microscopic calculation of the interaction cross section for stable and unstable nuclei 
based on the nonrelativistic nucleon-nucleon t-matrix 
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Fully quantal calculations of the total reaction cross sections ur and interaction cross sections 
(Ti, induced by stable and unstable He, Li, C and O isotopes on ^^C target at iSi^b ~ 0.8 and 1 
GeV/nucleon have been performed, for the first time, in the distorted wave impulse approximation 
(DWIA) using the microscopic complex optical potential and inelastic form factors given by the 
folding model. Realistic nuclear densities for the projectiles and ^^C target as well as the complex t- 
matrix parameterization of free nucleon-nucleon interaction by Franey and Love were used as inputs 
of the folding calculation. Our parameter-free folding + DWIA approach has been shown to give a 
very good account (within 1-2%) of the experimental ai measured at these energies for the stable, 
strongly bound isotopes. With the antisymmetrization of the dinuclear system properly taken into 
account, this microscopic approach is shown to be more accurate than the simple optical limit of 
Glauber model that was widely used to infer the nuclear radii from the measured ctj. Therefore, 
the results obtained for the nuclear radii of neutron-rich isotopes under study can be of interest for 
further nuclear structure studies. 

PACS numbers: 24.10.Eq, 24.10.Ht, 24.50.+g, 25.60.Bx, 25.60.Dz 



I. INTRODUCTION 

Since the 80s of the last century, the radioactive ion 
beams have been used intensively to measure the total 
reaction cross sections and interaction cross sections in- 
duced by unstable nuclei on stable targets (see a recent 
review in Ref. [1]) which serve as an important data bank 
for the determination of nuclear sizes. The discovery of 
exotic structures of unstable nuclei, such as neutron halos 
or neutron skins, are among the most fascinating results 
of this study. 

The theoretical tool used dominantly by now to ana- 
lyze the interaction cross sections measured at energies of 
several hundred MeV/nucleon is the Glauber model [2, 3] 
which is based on the eikonal approximation. This ap- 
proach provides a simple connection between the ground 
state densities of the two colliding nuclei and the total 
reaction cross section of the nucleus-nucleus system and 
has been used, in particular, to deduce the nuclear den- 
sity parameters for the neutron-rich halo nuclei [4]. 

In general, the total reaction cross section ctr, which 
measures the loss of flux from the elastic channel, must 
be calculated from the transmission coefficient Ti as 
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where k is the relative momentum (or wave number). 
The summation is carried over all partial waves / with T/ 
determined from the elastic S'-matrix as 

Ti = l-\Si\\ (2) 

In the standard optical model (OM), the quantal S'- 
matrix elements 5"/ are obtained from the solution of the 
Schrodinger equation for elastic nucleus-nucleus scatter- 
ing using a complex optical potential. At low energies, 
the eikonal approximation is less accurate and, instead 
of Glauber model, the OM should be used to calculate 
(TRfor a reliable comparison with the data. At energies 
approaching 1 GeV/nucleon region, there are very few 
elastic scattering data available and the choice of a re- 
alistic optical potential becomes technically difficult, es- 
pecially for unstable nuclei. Perhaps, that is the rea- 
son why different versions of Glauber model are widely 
used to calculate ctr at high energies. Depending on the 
structure model for the nuclear wave functions used in 
the calculation, those Glauber model calculations can be 
divided into two groups: the calculations using a simple 
optical limit of Glauber model (see Ref. [1] and refer- 
ences therein) and the more advanced approaches where 
the few-body correlation and/or breakup of a loosely- 
bound projectile into a core and valence (halo) nucleons 
are treated explicitly [3, 5, 6]. 

In the present work, we explore the applicability of the 
standard OM to calculate the total reaction cross section 
(1) induced by stable and unstable beams at high ener- 
gies using the microscopic optical potential predicted by 
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the folding model. The basic inputs of a folding calcula- 
tion are the densities of the two colliding nuclei and the 
effective nucleon-nucleon (NN) interaction [7]. At low 
energies, a realistic density dependent NN interaction [8] 
based on the M3Y interaction [9] has been successfully 
used to calculate the a-nucleus and nucleus- nucleus opti- 
cal potential [10]. This interaction fails, however, to pre- 
dict the shape of the a-nucleus optical potential as the 
bombarding energy increases to about 340 MeV/nucleon 
[11]. On the other hand, at incident energies approaching 
a few hundred MeV/nucleon the f-matrix parameteriza- 
tion of free NN interaction was often used in the folding 
analysis of proton-nucleus scattering [12, 13]. The use 
of the f-matrix interaction corresponds to the so-called 
impulse approximation (lA), where the medium modifi- 
cations of the NN interaction are neglected [14]. 

In the present folding calculation we adopt a local rep- 
resentation of the free NN f-matrix developed by Franey 
and Love [13] based on the experimental NN phase shifts. 
The folded optical potentials and inelastic form factors 
are used further in the distorted wave impulse approx- 
imation (DWIA) to calculate ctr and interaction cross 
section (Tj, induced by stable and unstable He, Li, C, 
and O isotopes on ^^C target at bombarding energies 
around 0.8 and 1 GeV/nucleon. Since relativistic effects 
are significant at high energies, the relativistic kinemat- 
ics are taken into account properly in both the folding 
and DWIA calculations. To clarify the adequacy and 
possible limitation of the present folding model, we also 
discuss the main approximations made in our approach 
and compare them with those usually assumed in the 
Glauber model. 

Given the realistic nuclear densities and validity of lA, 
the folding approach presented below in Sec. II is actually 
parameter-free and it is necessary to test first the relia- 
bility of the model by studying the known stable nuclei 
before going to study unstable nuclei. Such a procedure 
is discussed briefly in Sec. III. Then, ai measured for the 
neutron-rich He, Li, C, and O isotopes are compared with 
the results of calculation and the sensitivity of nuclear 
radii to the calculated ai is discussed. The discrepancy 
between o-p''^ and dj^^ found for some light halo nuclei 
is discussed in details to indicate possible effects caused 
by the dynamic few-body correlation. Conclusions are 
drawn in Sec. IV. 



II. FOLDING MODEL FOR THE COMPLEX 
NUCLEUS-NUCLEUS OPTICAL POTENTIAL 

The details of the latest double- folding formalism are 
given in Ref. [10] and we only recall briefly its main fea- 
tures. In general, the projectile-target interaction poten- 
tial can be evaluated as an energy dependent Hartree- 
Fock-type potential of the dinuclear system 

= ^ [< ij\vD\ij > + < ij\vEx\ji >] = Vd -|- Vex, 

(3) 



where the nuclear interaction y is a sum of effective NN 
interactions Vij between nucleon i in the projectile a and 
nucleon j in the target A. The antisymmetrization of 
the dinuclear system is done by taking into account the 
single-nucleon knock-on exchanges. 

The direct part of the potential is local (provided that 
the NN interaction itself is local), and can be written in 
terms of the one-body densities, 

Vb{E,R) = j pa{ra)pA{rA)vD{E,p,s)d^rc,d^rA, 

where s = va — Va + R- (4) 

The exchange part is, in general, nonlocal. However, an 
accurate local approximation can be obtained by treating 
the relative motion locally as a plane wave [15] 

Vex{E,R) = j Pa{ra,ra + s)pA{rA,rA - s) 

X VEx{E, p, s) exp (fracfrA. (5) 

Here Pa{ra) = Pa{ra, ra) and Pa{ra,ra + s) is the diago- 
nal and nondiagonal parts of the one-body density matrix 
for the projectile, and similarly for the target. K{E, R) 
is the local momentum of relative motion determined as 

K\E, R) = ^[£;c.m. - Re U{E, R) - Vc{R)], (6) 

p is the reduced mass, M = aA/ (a + A) with a and A the 
mass numbers of the projectile and target, respectively. 
Here, U{E,R) = Vu{E,R) + Vex{E,R) and Vc{R) are 
the total nuclear and Coulomb potentials, respectively. 
More details on the calculation of the direct and exchange 
potentials (4)- (5) can be found in Refs. [10, 16]. The fold- 
ing inputs for mass numbers and incident energies were 
taken as given by the relativistically corrected kinematics 
[17]. 

To calculate consistently both the optical potential and 
inelastic form factor one needs to take into account ex- 
plicitly the multipole decomposition of the nuclear den- 
sity that enters the folding calculation [10] 

PjM^j'M'ir) = ^ < JMXp\J'M' > CxPx{r)[i^Yx4f)]*, 

(7) 

where JM and J'M' are the nuclear spin and its pro- 
jection in the initial and final states, respectively, and 
P\{r) is the nuclear transition density for the correspond- 
ing 2'^-pole excitation. In the present work, we adopt the 
collective-model Bohr-Mottelson prescription [18] to con- 
struct the nuclear transition density for a given excitation 
in the ^^C target as 

dpo{r) 

^aW = (8) 

Here po{r) is the total ground state (g.s.) density and 6x 
is the deformation length of the 2^-pole excitation in the 
i^C target. 
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Impulse approximation and the f-matrix interaction 

If the total spin and isospin are zero for one of the 
two colliding nuclei ("'^C in our case) only the spin- and 
isospin-independent components of the central NN forces 
are necessary for the folding calculation. We discuss now 
the choice of wd(ex) {E, p, s) for the two bombarding en- 
ergies of 0.8 and 1 GeV/nuclcon. At these high energies, 
one can adopt the lA which reduces the effective NN in- 
teraction approximately to that between the two nucleons 
in vacuum [14]. Consequently, the microscopic optical 
potential and inelastic form factors can be obtained by 
folding the g.s. and transition densities of the two collid- 
ing nuclei with an appropriate t-matrix parameterization 
of the free NN interaction. 

In the present work, we have chosen the nonrelativis- 
tic t-matrix interaction which was developed by Praney 
and Love [13] based on experimental NN phase shifts 
at bombarding energies of 0.8 and 1 GeV. The spin- and 
isospin-independent direct (wd) and exchange (uex) parts 
of the central NN interaction are then determined from 
the singlet- and triplet-even (SE,TE) and odd (SO,TO) 
components of the local t-matrix interaction (see Table I 
of Ref. [13]) as 
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[3iTE (S) + 3tsE (S) ± 9tTO (s) ± 3tso (s) 



(9) 

Here ka and kA are the energy-dependent kinematic mod- 
ification factors of the t-matrix transformation [19] from 
the NN frame to the Na and NA frames, respectively, ka 
and kA were evaluated using Eq. (19) of Ref. [12]. The 
explicit, complex strength of the finite-range central t- 
matrix interaction (9) is given in terms of four Yukawas 
[13]. Since the medium modifications of the NN interac- 
tion are neglected in the lA [14] , the t-matrix interaction 
(9) does not depend on the nuclear density. 



Main steps in the calculation of or 

With properly chosen g.s. densities for the two collid- 
ing nuclei, the elastic scattering cross section and (Tr are 
obtained straightforwardly in the OM calculation using 
the microscopic optical potential (4)- (6). We recall that 
the interaction cross section ci is actually the sum of all 
particle removal cross sections from the projectile [1] and 
accounts, therefore, for all processes when the neutron 
and/or proton number in the projectile is changed. As a 
result, (Ti must be smaller than the total reaction cross 
section ctr wliic;li includes also the cross section of in- 
elastic scattering to excited states in both the target and 
projectile as well as cross section of nucleon removal from 
the target. At energies of several hundred MeV/nucleon, 
the difference between ur and ai was found to be a few 
percent [3, 20, 21] and was usually neglected to allow a 
direct comparison of the calculated ctr with the measured 
(7i. Since the experimental uncertainty in the measured 



(7i is very small at the considered energies (around 1% for 
stable projectiles like ^He, ^^C, and ^^O [1]) neglecting 
the difference between ctr and ui might be too rough an 
approximation in comparing the calculated ctr with the 
measured ai and testing nuclear radius at the accuracy 
level of ±0.05 fm or less [1, 22]. In the present work, 
we try to estimate ai as accurately as possible by sub- 
tracting from the calculated <tr the total cross section of 
the main inelastic scattering channels. Namely, we have 
calculated in DWIA, using the complex folded optical 
potential and inelastic form factors, the integrated cross 
sections 172+ ^-nd (T3- of inelastic scattering to the first 
excited 2+ and 3^ states of ^^C target at 4.44 and 9.64 
MeV, respectively. These states are known to have the 
largest cross sections in the inelastic proton and heavy 
ion scattering on ^^C at different energies. The deforma- 
tion lengths used to construct transition densities (8) for 
the folding calculation were chosen so that the electric 
transition rates measured for these states are reproduced 
with the proton transition density as 



B{EX t) = 



,A+2 



dr 



(10) 



Using a realistic Fermi distribution for the g.s. den- 
sity of ^^C (see next Section) to generate the transi- 
tion densities, we obtain 62 ~ 1.54 fm and ^3 « 2.11 
fm which reproduce the experimental transition rates 
B{E2 t) « 41 e^fm" [23] and B{E[i t) ~ 750 e^fm^ [24], 
respectively, via Eq. (10). Since inelastic scattering to 
excited states of the unstable projectile is suppressed by 
a much faster breakup process, ai can be approximately 
obtained as 



CTI = CTR — CTlnel 

» (Tr — C72+ — (73- 



(11) 



All the OM and DWIA calculations were made using 
the code ECIS97 [25] with the relativistic kinematics 
properly taken into account. At the energies around 1 
GeV/nucleon the summation (1) is usually carried over 
up to 800 - 1000 partial waves to reach the full conver- 
gence of the S'-matrix series for the considered nucleus- 
nucleus systems. 



Adequacy and limitation of the folding approach 

Since the measured txi have been analyzed extensively 
by different versions of Glauber model and its optical 
limit (OL) is sometimes referred to as the folding model 
[6, 26], we find it necessary to highlight the distinctive 
features of the present folding approach in comparison 
with the OL of Glauber model before going to discuss 
the results of calculation. 

On the level of the nucleus-nucleus optical potential 
(OP), the present double- folding approach evaluates OP 
using fully finite-range NN interaction and taking into 
account the exchange effects accurately via the Fock term 
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FIG. 1: Radial shape of the direct Vr> and exchange Vex parts of the total optical potential U for ®He+^^C system at 790 
MeV/nucleon. The real and imaginary parts of U are shown in panels a and b, respectively. 



in Eq. (3). Therefore, individual micleons are allowed to 
scatter after the collision into unoccupied single-particle 
states only. Sometimes, one discusses these effects as the 
exchange NN correlation. An appropriate treatment of 
the exchange NN correlation is indispensable not only in 
the folding calculation of OP and inelastic form factor, 
but also in the Hartree-Fock (HF) calculations of nuclear 
matter [27] and of the finite nuclei [28]. 

To obtain from the double-folding model presented 
above the simple expression of nucleus-nucleus OP used 
in the OL of Glauber model one needs to make a "double- 
zero" approximation which reduces the complex finite- 
range t-matrix interaction (9) to a zero-range (purely 
imaginary) NN scattering amplitude at zero NN angle 
iNN(^ = 0°)(5(s) that can be further expressed through 
the total NN cross section ctnn, using the optical the- 
orem. As a result, one needs to evaluate in the OL of 
Glauber model only a simple folding integral over local 
densities of the two colliding nuclei [6] 

UiR) ^ Vol{R) = !^ y p,{R)pA{R-rA)d^rA. 

(12) 

The prescription (12) is also known as the tpp approx- 
imation [29] which neglects the off-shell part of the t- 
matrix. Besides the inaccuracy caused by the use of zero- 
range approximation [30], the zero-angle approximation 
takes into account only the on-shell t-matrix at zero mo- 
mentum transfer (sec Eq. (3) in Rcf. [12]). Since the an- 
tisymmetrization of iNN requires an accurate estimation 



of the NN knock-on exchange term which is strongest 
at large momentum transfers (q > 6 fm^^ at energies 
around 0.8 GcV [12, 13]), the zero-angle approximation 
could strongly reduce the strength of the exchange term. 
A question remains, therefore, whether the NN antisym- 
metry is properly taken into account when one uses the 
empirical o-nn in the Glauber folding integral (12). A 
similar aspect has been raised by Brandan et al. [31] who 
foimd that an overestimated absorption in the nucleus- 
nucleus system (by the tpp model) is due to the effects 
of Pauli principle. To illustrate the importance of the 
knock-on exchange term, we have plotted in Fig. 1 the 
direct and exchange components of the microscopic OP 
for ^He+^^C system at 790 MeV/nucleon predicted by 
our double-folding approach using realistic g.s. densities 
(see next Sect.) of the two colliding nuclei. One can sec 
that the exchange term of the real OP is repulsive and 
much stronger than the (attractive) direct term, which 
makes the total real OP repulsive at all internuclear dis- 
tances (sec panel a of Fig. 1). The exchange term of the 
imaginary OP is also repulsive but its relative strength 
is much weaker compared to that of the real OP, and the 
total imaginary OP remains attractive or absorptive at 
all distances. As a result, the direct part of the imag- 
inary OP is about 10% more absorptive than the total 
imaginary OP (see panel b of Fig. 1). The total reac- 
tion cross section predicted by the complex OP shown in 
Fig. 1 is (Tr « 727 mb. This value increases to ctr w 750 
mb when the exchange potential Vex is omitted in the 
OM calculation. Consequently, the relative contribution 
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by the exchange term in ctr is about 3%. This difference 
is not small because it can lead to a difference of up to 
7% in the extracted nuclear RMS radii. Due to an over- 
whelming contribution by the exchange part of the real 
OP, the exchange potential affects the calculated elastic 
scattering cross section (sec; Fig. 2) much more substan- 
tially compared to ctr., which is determined mainly by 
the imaginary OP. 

We will show below a slight (but rather systematic) 
difference in ctr values obtained in our approach and 
the OL of Glauber model that might be due to the ex- 
change effect. Wc note hirthcr that the elastic S'-matrix 
is obtained in our approach rigorously from the quantal 
solution of the Schrodinger equation for elastic scatter- 
ing wave, while the clastic S'-matrix used in the Glauber 
model is given by the eikonal approximation which ne- 
glects the second-derivative term of the same Schrodinger 
equation. 

A common feature of the present folding approach and 
the OL of Glauber model is the use of single-particle nu- 
clear densities of the projectile and target as input for the 
calculation, leaving out all few-body correlations to the 
structure model used to construct the density. This sim- 
ple ansatz has been referred to as "static density approxi- 
mation" [5, 6] which does not take into account explicitly 
the dynamic few-body correlation between the core and 
valence nucleons in a loosely-bound projectile while it 
collides with the target. In the Glauber model, this type 
of few-body correlation can be treated explicitly [3, 5, 6] 
using simple assumptions for the wave functions of the 
core and valence nucleons as well as that of their relative 
motion. For unstable nuclei with a well-extended halo 
structure, like ^^Li or ^He, such an explicit treatment of 
the dynamic few-body correlation leads consistently to a 
smaller ctr , i.e., to a larger nuclear radius compared to 
that given by the OL of Glauber model [3, 5, 6]. On 
the level of the HF-type folding calculation (3), an ex- 
plicit treatment of the core and valence nucleons would 
result in a much more complicated triple-folding formal- 
ism which involves the antisymmetrization not only be- 
tween the projectile nucleons and those of the target, but 
also between the nucleons of the core and the valence nu- 
cleons. Such an approach would clearly end up with a 
nonlocal OP which will not be easily used with the ex- 
isting direct reaction codes. The lack of an appropriate 
treatment of the dynamic few-body correlations remains, 
therefore, the main limitation of the present folding ap- 
proach in the calculation of the OP for systems involving 
unstable nuclei with halo-type structure. 

Note that an effective way of taking into account the 
loose binding between the core and valence nucleons is to 
add a higher-order contribution from breakup (dynamic 
polarization potential) to the first-order folded potential 
[21, 32] or simply to rcnormalize the folded potential to 
fit the data. However, validity of the lA implies that 
higher-order multiple scattering or contribution from the 
dynamic polarization potential are negligible, and the 
folded OP and inelastic form factor based on the f-matrix 



interaction (9) should be used in the calculations with- 
out any further renormalization. Therefore, we will dis- 
cuss below only results obtained with the unrenormalized 
folded potentials, keeping in mind possible effects due to 
the few-body correlation. 

III. RESULTS AND DISCUSSION 

Results for stable (A'^ = Z) isotopes 

An important step in any experimental or theoreti- 
cal reaction study with unstable beams is to gauge the 
method or model by the results obtained with stable 
beams. Therefore, wc have considered first the available 
data of CTi induced by stable ''He, ^Li, ^^C, and ^^O beams 
on ^^C target [1]. These {N = Z) nuclei are strongly 
bound, and the RMS radius of the (point) proton distri- 
bution inferred from the elastic electron scattering data 
[35] can be adopted as the "experimental" nuclear radius 
if the proton and neutron densities are assumed to be 
the same. To show the sensitivity of the calculated ai to 
the nuclear radius, we present in Table I results obtained 
with different choices for the projectile density in each 
case. We use for the g.s. density of ^^C target a realistic 
Fermi (FM) distribution [16] 

Po(r-) = Pol [I + exp((r - c)/a)], (13) 

where po = 0.194 fm^^^ ^ ^ 2.214 and a = 0.425 fm were 
chosen to reproduce the shape of shell model density and 
experimental radius of 2.33 fm for ^^C. 

^He is a unique case where a simple harmonic oscil- 
lator (HO) model can reproduce quite well its ground 
state density. If one chooses the HO parameter to give 
< >i/2= 1.461 fm (close to the experimental radius 
of 1.47 ± 0.02 fm) then one obtains the Gaussian form 
adopted in Ref. [7] for a-density. This choice of '*He 
density has been shown in the folding analysis of elas- 
tic a-nucleus scattering [16] to be the most realistic. By 
comparing the calculated cti with the data, we find that 
this same choice of ''He density gives the best agreement 
between ctj*''^ and ctj^^. Similar situation was found for 
^^C and '^^O isotopes, where the best agreement with the 
data is given by the densities which reproduce the exper- 
imental nuclear radii. Beside a simple Fermi distribution 
[16], microscopic g.s. densities given by the Hartree-Fock- 
Bogoliubov (HFB) calculation that takes into account 
the continuum [33] were also used. The agreement with 
the data for ^^C and ^^O given by the HFB densities 
is around 2%, quite satisfactory for a fully microscopic 
structure model. We have further used sp-shcll HO wave 
functions to construct the g.s. densities of ^Li, '^^C, and 
^^O. For ^^C and ^^O, the best agreement with the ai 
data is again reached when the HO parameter is tuned 
to reproduce the experimental radii. 

The agreement is slightly worse for ^Li compared to 
^He, ^^C, and ^^O cases if ''Li density distribution repro- 
duces the experimental radius. We have first used ^Li 
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TABLE I: The total reaction cross section ctr and interaction cross section ai calculated for stable ''He, ^''^Li, ^^C, and 
nuclei in comparison with cr^^^ taken from the data compilation in Ref. [1]. Aai = IctP'"^ — crj''''|/crj'^^. 



Nucleus 


Energy 
(MeV/nucleon) 


Density Model 


(fm) 


Reference 


< >J/^ 
(fm)"' 


calc 
tjR 

(mb) 


calc 

(mb) 


cxp 

(mb) 


Acri 

(%) 


*He 


790 


HO 


1.461 


[7] 


1.47 ±0.02 " 


513 


504 


503 ±5 


0.2 






HO 


1.550 


[16] 




523 


515 




2.4 






HO 


1.720 


[36] 




543 


536 




6.6 


"Li 


790 


IPM 


2.401 


[34] 


2.43 ± 0.02 " 


722 


717 


688 ± 10 


4.2 






HO 


2.401 


This work 




723 


718 




4.4 






HO 


2.320 


This work 




709 


703 




2.2 


'Li 


790 


IPM 


2.367 


[34] 


2.33 ± 0.02 " 


746 


741 


736 ±6 


0.7 






HO 


2.334 


This work 




744 


739 




0.4 




950 


FM 


2.332 


[16] 


2.33 ±0.02 " 


854 


844 


853 ±6 


1.1 






HO 


2.332 


[16] 




853 


843 




1.1 






HFB 


2.446 


This work 




881 


872 




2.2 




970 


FM 


2.618 


[16] 


2.61 ±0.01 " 


992 


981 


982 ±6 


0.1 






HO 


2.612 


[16] 




988 


978 




0.4 






HFB 


2.674 


This work 




1006 


997 




1.4 



" RMS radius of the proton density given by the experimental charge density [35] unfolded with the finite size of proton. 
^ Nuclear RMS radius deduced from the Glauber model analysis of the same ai data in the OL approximation [1] . 



density given by the independent particle model (IPM) 
developed by Satchler [7, 34] which generates realistic 
wave function for each single-particle orbital using a 
Woods-Saxon (WS) potential for the bound state prob- 
lem. The IPM density gives < >i/2« 2.40 fm for 
®Li, rather close to the experimental radius of 2.43 ±0.02 
fm inferred from (e,e) data [35]. The HO density gives 
the same cti as that given by the IPM density if the HO 
parameter is chosen to give the same radius of 2.40 fm. 
These two versions of ^Li density overestimate the ai 
data by about 4%. If the HO parameter is chosen to 
give < 2.32 fm then the agreement with the 

ai data improves to around 2%. This result indicates 
that our folding + DWIA analysis slightly overestimates 
the absorption in ^Li+^^C system. Since ^Li is a loosely 
bound a + d system, this few percent discrepancy with 
the fTi data might well be due to the dynamic correlation 
between the a-corc and dcutcron cluster in ^Li during 
the collision which is not taken into account by our ap- 
proach. Note that a few-body Glauber calculation [6] 
(which takes into account explicitly the dynamic correla- 
tion between a and d) ends up, however, with about the 
same discrepancy (see Fig. 4 in Ref. [6]). ^Li remains, 
therefore, an interesting case for the reaction models to 
improve their ingredients. For "^lA, the IPM density [7] 
gives < >p^^» 2.28 fm (close to the experimental value 
of 2.27 ± 0.01 fm [35]) and < >y^« 2.43 fm which 
make the matter radius < >^/^f« 2.37 fm. As a result. 
(Ti calculated with the IPM density for ''Li agrees with 
the data within less than 1%. In the HO model for ''Li 
density, we have chosen the HO parameter for protons to 
reproduce the experimental radius of 2.27 fm and that 
for neutrons adjusted by the best agreement with the ai 
data. The best-fit < >^/^radius then becomes around 
2.33 fm. 



We conclude from these results that the present fold- 
ing + DWIA approach and local t-matrix interaction by 
Praney and Love [13] are quite suitable for the descrip- 
tion of the nucleus-nucleus interaction cross section at 
energies around 1 GeV/nucleon, with the prediction ac- 
curacy as fine as 1-2% for the stable and strongly bound 
nuclei. 



Results for neutron-rich isotopes 

Our results for neutron-rich He, Li, C, and O isotopes 

arc presented in Table II. Since ^He beams are now 
available with quite a good resolution, this nucleus is 
among the most studied unstable nuclei. In the present 
work we have tested 3 different choices for ^He density 
in the calculation of cti. The microscopic ^He density ob- 
tained in a HF calculation [30] has a rather small radius 
< j,2 >i/2^ 2.20 fm and the calculated ui underestimates 
the data by about 5%. A larger radius of 2.53 fm is given 
by the density obtained in a consistent three-body for- 
malism [5] and the corresponding ai agrees better with 
the data. Given an accurate '^Li density obtained in the 
IPM [7] as shown above and the fact that ^He can be 
produced by a proton-pickup reaction on 'Li, we have 
constructed the g.s. density of ^He in the IPM (with the 
recoil effect properly taken into account [34]) using the 
following WS parameters for the single-particle states: 
ro = 1.25 fm, a = 0.65 fm for the neutrons and pro- 
tons which are bound by Sn = 25 MeV and Sp = 23 
MeV, respectively; ro = 1.35 fm, a = 0.65 fm for the p| 
halo neutrons which are bound by Sn = 1-86 MeV. The 
WS depth is adjusted in each case to reproduce the bind- 
ing energy. The obtained IPM density gives the proton, 
neutron and total nuclear radii of ^He as 1.755, 2.746 
and 2.460 fm, respectively. This choice of ^He density 
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TABLE II: The same as Table I but for neutron-rich He, Li, C, and O isotopes. Note that < >caic given by the HO densities 
should have about the same uncertainties as those deduced for < >llp by the OL of Glauber model. 
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" Nuclear RMS radius deduced from the Glauber model analysis of the ai data in the OL approximation [1] . 
^ Nuclear RMS radius deduced from the Glauber model analysis of elastic ^He+p scattering data at 0.7 GeV/nucleon [38]. 

° (Ti data taken from Ref . [43] . 



also gives the best agreement with the ai data. We note 
that a Glauber model analysis of the elastic ^He+p scat- 
tering at 0.7 GeV/nucleon [38], which takes into account 
higher-order multiple scattering effects, gives a best- fit 
< >i/^« 2.45 fm for ^He, very close to our result. 
Since elastic ^He+^'^^C scattering has recently been mea- 
sured at lower energies [40] , we found it interesting to plot 
the 3 densities and elastic ^He+^^C scattering cross sec- 
tions at 790 MeV/nucleon predicted by the corresponding 
complex folded OP (the radial shape of the OP obtained 
with the IPM density for ^He is shown in Fig. 1). As can 
be seen from Fig. 2, the IPM density has the neutron- 
halo tail very close to that of the density calculated in the 
three-body model [5] and they both give a good descrip- 



tion of (Tj. The predicted elastic cross section is strongly 
forward peaked and the difference in densities begins to 
show up after the first diffractive maximum. Such a mea- 
surement should be feasible at the facilities used for elas- 
tic ^He+p scattering at 0.7 GeV/nucleon [38] and would 
be very helpful in testing finer details of ^Hc density. As 
already discussed in previous Sect., the exchange part of 
the microscopic OP affects the elastic cross section very 
strongly (see dotted curve in panel h of Fig. 2) and the 
elastic ^He+^^C scattering measurement would be also a 
very suitable probe of the exchange effects in this system. 

Since ^ He is a loosely bound halo nucleus with a well es- 
tablished three-body a+n+n structure, the dynamic cor- 
relation between the a-core and dineutron is expected to 
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FIG. 2: Three versions of ^He g.s. density used in the folding calculation (panel a) and elastic ®He+^^C scattering cross sections 
at 790 McV/nuclcon obtained with the corresponding complex folded optical potentials (panel 6). The dotted curve in panel 
b is obtained without the exchange part of the OP. 



be important during the collision. Our folding + DWIA 
approach using three-body density for ^He (version FC 
[5]) gives ai « 733 mb compared to about 720 mb given 
by the few-body calculation by Tostevin et al. (see Fig. 4 
in Ref. [6]) based on the same three-body wave func- 
tion for ^He. The difference in the calculated cti leads 
to an increase of about 2-3% in the < >^/^ value. It 
is likely that such a difference is, in part, due to the 
dynamic correlation between the a-core and dineutron 
which was not considered in our folding -|- DWIA ap- 
proach. For ®He nucleus, the OL of Glauber analysis of 
Cti data [1] and the multiple-scattering Glauber analy- 
sis of elastic ®He+p data at 0.7 GeV/nucleon [38] give 
< >^/^around 2.52 and 2.53 fm, respectively. By us- 
ing the microscopic ^He density obtained in a four-body 
(COSMA) model [37], which gives < >^/'^= 2.526 fm, 
our folding -|- DWIA approach reproduces the measured 
Cti data within less than 1%. Note that a (multiple scat- 
tering) Glauber model analysis of the elastic ^'*He+p 
scattering at 0.7 GeV/nucleon which takes into account 
the dynamic few-body correlation explicitly was done by 
Al-Khalili and Tostevin [39] , and they have obtained the 
best-fit nuclear radii of about 2.5 and 2.6 fm for ^He and 
^He, respectively, around 2% larger than our results. 



Parameters of HO densities deduced from ai data 

Although the HO model is a very simple approach, 
the HO densities were shown above to be useful in test- 
ing the nuclear radii for stable (A'" = Z) nuclei. More- 
over, the HO-type densities (with the appropriately cho- 
sen HO lengths) for the sd-shell nuclei have been suc- 
cessfully used in the analysis of (e,e) data, measurements 
of isotope shift and muonic atoms [1]. Therefore, it is 
not unreasonable to use simple HO parameterization for 
the g.s. densities of neutron-rich nuclei to estimate the 
nuclear radii, based on our folding -|- DWIA analysis of 
Cti data. For a N ^ Z nucleus, one needs to generate 
proton and neutron densities separately as 

Mr) = ^ (l + P.g + ^.g) exp (-0 , (14) 

where r = n or p, parameters Pr and Dr are determined 
from the nucleon occupation of the p- and d harmonic- 
oscillator shells, respectively. 

To generate the g.s. densities of ^'^Li isotopes, we have 
assumed the proton density of these nuclei to be approx- 
imately that of ''Li and the neutron HO length 6„ is ad- 
justed in each case to reproduce the measured cti (see 
Tables II and HI). While the obtained < >^/Hot «Li 
is rather close to that given by the OL of Glauber model 
[1], results obtained for ''Li arc different and we could 
reproduce the cti data only if the neutron HO length is 
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TABLE III: The HO-density parameters (14) for neutron-rich Li, C, and O isotopes. 



Nucleus 


Pn 


Pp 




Dp 


bn 


bp 


< >M^ 


< >T 














(fm) 


(fm) 


(fm) 


(fm) 


(fm) 


'^Li 


2/3 


1/3 


0.0 


0.0 


1.684 


1.6766 


2.382 


2.270 


2.334 


8Li 


1.0 


1/3 


0.0 


0.0 


1.6770 


1.6776 


2.430 


2.270 


2.371 


^Li 


4/3 


1/3 


0.0 


0.0 


1.6470 


1.6766 


2.424 


2.270 


2.374 


13 c 


5/3 


4/3 


0.0 


0.0 


1.6058 


1.5722 


2.389 


2.314 


2.355 


"C 


2.0 


4/3 


0.0 


0.0 


1.6226 


1.5762 


2.434 


2.320 


2.386 




2.0 


4/3 


2/15 


0.0 


1.6630 


1.5898 


2.570 


2.340 


2.481 




2.0 


4/3 


4/15 


0.0 


1.8512 


1.7128 


2.927 


2.521 


2.782 




2.0 


4/3 


2/5 


0.0 


1.8552 


1.7128 


2.986 


2.521 


2.831 




2.0 


4/3 


8/15 


0.0 


1.8752 


1.7297 


3.062 


2.546 


2.900 




2.0 


4/3 


2/3 


0.0 


2.1252 


1.7297 


3.512 


2.546 


3.238 


20 c 


2.0 


4/3 


4/5 


0.0 


1.9462 


1.7467 


3.248 


2.571 


3.061 


^'O 


2.0 


2.0 


2/15 


0.0 


1.7775 


1.7232 


2.747 


2.585 


2.672 


18q 


2.0 


2.0 


4/15 


0.0 


1.7601 


1.7935 


2.783 


2.690 


2.742 


19 Q 


2.0 


2.0 


2/5 


0.0 


1.7601 


1.7935 


2.833 


2.690 


2.774 


20q 


2.0 


2.0 


8/15 


0.0 


1.7401 


1.8005 


2.842 


2.701 


2.786 


21 Q 


2.0 


2.0 


2/3 


0.0 


1.7401 


1.8005 


2.876 


2.701 


2.811 


22q 


2.0 


2.0 


4/5 


0.0 


1.8498 


1.8081 


3.087 


2.712 


2.956 


23 Q 


2.0 


2.0 


14/15 


0.0 


2.1118 


1.8081 


3.555 


2.712 


3.286 


24q 


2.0 


2.0 


16/15 


0.0 


2.0758 


1.8261 


3.520 


2.739 


3.280 



chosen to give < >caic~ "^'^"^ ^™ about 2% larger 
than that given by the OL of Glauber model. For the 
halo nucleus ^^Li, a ^Li core + two-neutron halo model 
was used to generate its density. Namely, we have used 
HO density of ^Li that reproduces the measured <j\ for 
^Li and a Gaussian tail for the two-neutron halo den- 
sity. To reach a best agreement between taken from 
Ref. [43] and o-p'^, the Gaussian range was chosen to give 

< >ciic~ "^-^"^ ^™ which is about 0.1 fm larger than 
that given by the OL of Glauber model [1]. A micro- 
scopic density for '^^Li obtained in the HF calculation [30] 
(which gives < >i/2= 2.868 fm) has also been used 
in our folding analysis. The agreement with the data 
becomes much worse in this case (see Table II) and we 
conclude that the radius given by the HF density is some- 
what too small. To show the sensitivity of our analysis to 
the nuclear radius, we have plotted in Fig. 3 a\ predicted 
by three versions of ^^Li density with the Gaussian range 
of the 2n-halo adjusted to give < >i/2= 3.15^ 3.23 
and 3.30 fm, respectively, compared to = 1060 ± 10 
mb [43]. It is easily to infer from Fig. 3 an empirical 
RMS radius of 3.23 ±0.05 fm for "Li. Note that txi mea- 
surement for ^^Li-h^^C system at 790 MeV/nucleon has 
been reported in several works with = 1040 ± 60 
[42], 1047 ± 40 [41] and 1060 ± 10 mb [43]. If we adjust 
Gaussian range of the 2n-halo in ^^Li density to repro- 
duce these a^^^ values, the corresponding < >^/^radii 
of -'^'^Li are 3.13, 3.15 and 3.23 fm, respectively. Since 
(Ti data obtained in Ref. [43] has a much better statis- 
tics and less uncertainty, we have adopted < >i/2= 
3.23 ± 0.05 fm as the most realistic RMS radius of ^^Li 
given by our folding -I- DWIA analysis. 

The total reaction cross section for ^^Li-|-^^C system 



at 790 McV/nucleon has been studied earlier in the few- 
body Glauber formalism by Al-Khalili et al. [5], where 
< >^/^radius for ^^Li was shown to increase from 
3.05 fm (in the OL) to around 3.53 fm when the dy- 
namic correlation between ^Li-core and 2n-halo during 
the collision is treated explicitly. This is about 9% larger 
than < >^ /^radius obtained in our folding -|- DWIA 
approach based on the same cti data. Although vari- 
ous structure calculations for ^^Li give its RMS radius 
around 3.1-3.2 fm (see Refs. [1, 4] and references therein), 
a very recent coupled-channel three-body model for "'^-'^Li 
by Ikeda et al. [44, 45] shows that its RMS radius is 
ranging from 3.33 to 3.85 fm if the 2n-halo wave function 
consists of 21 to 39% mixture from (s^)^ state, respec- 
tively. A comparison of the calculated Coulomb breakup 
cross section with the data [45] suggests that this s-wave 
mixture is around 20-30 %. Thus, the nuclear radius of 
^^Li must be larger than that accepted sofar [1, 4] and be 
aromid 3.3 3.5 fm, closer to the result of the few-body 
calculation [5] and the upper limit of RMS radius given 
by our folding -|- DWIA analysis. 



For most of neutron-rich C and O isotopes considered 
here, we have first fixed the proton HO lengths bp to 

reproduce the proton < >p^^radii predicted by the 
microscopic IPM and HFB densities (as described below). 
The neutron HO lengths 6„ are then adjusted to the best 
agreement with ai data, and the obtained HO parameters 
are summarized in Table HI. 
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FIG. 3: af' obtained with three versions of Li g.s. density, 
where Gaussian range of the 2'fi-halo was adjusted to give 
< >^^^=3.15, 3.23 and 3.30 fm for '^'^Li, in comparison 
with ffj^'"' = 1060 ± 10 mb [43]. 



Microscopic HFB densities 

Before discussing the results obtained for the neutron- 
rich C and O isotopes, wc give here a brief description of 
the microscopic HFB approach used to calculate the g.s. 
densities of even C and O isotopes. More details about 
this approach can be found in Ref. [33]. 

We solve the HFB equations in coordinate represen- 
tation and in spherical symmetry with the inclusion of 
continuum states for neutron-rich nuclei. As the neu- 
tron Fermi energies of these nuclei are typically quite 
close to zero, pairing correlations can easily scatter pairs 
of neutrons from the bound states towards continuum 
states. For this reason, the inclusion and the treatment 
of continuum states in the calculation are very impor- 
tant. In our calculation the continuum is treated exactly, 
i.e., with the correct boundary conditions for continuum 
wave functions and by taking into account the widths of 
the resonances. Resonant states are localized by study- 
ing the behavior of the phase shifts with respect to the 
quasi-particle energy for each partial wave 

The calculations were done with the Skyrme interac- 
tion SLy4 for the mean field channel and with the follow- 



ing zero-range density-dependent interaction 



V = Vu 



1 



Pi) 



(5(ri - r2) 



(15) 



for the pairing channel. In Eq. (15), po is the saturation 
density and 7 is chosen equal to 1. We have adapted 
the prescription of Refs. [46, 47] to finite nuclei in or- 
der to fix Vq together with the quasi-particle energy cut- 
off. This prescription, requiring that the free neutron- 
neutron scattering length has to be reproduced in the 
truncated space, allows us to deduce a relation between 
the parameter Vq and the quasiparticle energy cutoff. 



Nuclear radii of Carbon and Oxygen isotopes 

The CTi data for neutron-rich C and O isotopes are com- 
pared in Table 11 with ai predicted by different choices 
of nuclear densities. We have tested first the 1PM den- 
sity for ^''C [34] based on the single-particle spectroscopic 
factors obtained in the shell model by Cohen and Kurath 
[48]. This IPM density gives < >^^^Ki 2.39 fm for 
and the predicted ai agrees with the data within less than 
2%. We have further made 1PM calculation for ^^C based 
on the same single-particle configurations, with the WS 
parameters for sp-shells appropriately corrected for the 
recoil effects and experimental nucleon separation ener- 
gies Sn,p of i^C. This 1PM density gives < >i/2- 2.42 
fm for ^"'C and the predicted cri also agrees with the data 
within 2%. The HO densities were also parameterized 
for 13,i4q ^j^j^ ^jig proton HO lengths bp chosen to repro- 
duce < >y^values predicted by the 1PM. The best-fit 
neutron HO lengths 6„ result in < >i/2= 2.36 and 
2.39 fm for ^'^C and ^''C, respectively. These values agree 
fairly with those given by the 1PM densities. The mi- 
croscopic HFB density gives for ^^C a significantly larger 
< >^/^radius of 2.59 fm and the calculated cti over- 
estimates the data by nearly 7%. Note that the OL of 
Glauber model gives smaller radius of 2.28 and 2.30 fm 
for ^'^C and ^"'C, respectively, based on the same cti data 
[1]. This means that the absorption given by the OL of 
Glauber model is indeed stronger than that given by our 
approach, as expected from discussion in Sect. II. 

For the neutron-rich even 16-20Q isotopes, the HFB 
densities give a remarkably better agreement with the 
data and it is, therefore, reasonable to fix the proton HO 
lengths of the HO densities for each of i5-20(~; isotopes to 
reproduce < >y^radius predicted by the HFB calcu- 
lation for the nearest even neighbor. The best-fit neutron 
HO lengths result in the nuclear radii quite close to those 
given by the HFB densities (see Tables II and III). We 
emphasize that the nuclear radii given by our analysis, 
using the HO densities for C isotopes, are about 0.1 fm 
larger than those deduced from the OL of Glauber model 
[1]. Given a high sensitivity of ai data to the nuclear size, 
a difference of 0.1 fm is not negligible. 

To illustrate the mass dependence of the nuclear ra- 
dius, we have plotted in Fig. 4a the RMS radii given 
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FIG. 4: Mass dependence of the nuclear RMS radius for Carbon (panel a) and Oxygen (panel b) isotopes given by the two 
choices (HFB and HO) of the g.s. densities compared to that deduced from the Glauber model analysis in the OL approximation 
[1]. The lines represent roA^^^ dependence with ro deduced from the experimental radii of ^^C and ^''O given in Table I. 



by the two sets (HFB and HO) of the g.s. densities 
for C isotopes together with those deduced from the OL 
of Glauber model based on the same cri data [1]. One 
can sec that our result follows closely the trend estab- 
lished by the OL of Glauber model, although the abso- 
lute < >^/^radii obtained with the HO densities are 
in most cases larger than those deduced from the OL of 
Glauber model. With the exception of the ^^C case, the 
radii of even C isotopes given by the microscopic HFB 
densities agree reasonably well with the empirical HO re- 
sults. We have also plotted in Fig. 4 the lines representing 
roA^/^ dependence with tq deduced from the experimen- 
tal radii of ^^C and ^^O given in Table I. One can see that 
the behavior of nuclear radius in C isotopes is quite dif- 
ferent from the roA^/'^ law. While < >^/^radii found 
for i2-i5(^ agj-QQ fairly with the vqA^^^ law, those ob- 
tained for 16-20 Q significantly higher. In particular, 
a jump in the < >^/^value was found in ^^C compared 
to those found for -'^^"^^C. This result seems to support 
the existence of a neutron halo in ^^C as suggested from 
the ctr measurement for this isotope at 85 MeV/nucleon 
[49]. We have further obtained a nuclear radius of 3.24 
fm for ^^C which is significantly larger than that found 
for ^°C. This result might also indicate to a neutron halo 
in this odd G isotope. 

Situation is a bit different for O isotopes, where the 
best-fit < >^/^radii follow roughly the roA^/^ law up 
to ^^O. For the stable ^''^i^q jgotopcs, the IPM densities 
[34] provide a very good description of the ai data (within 



1-2%). The best-fit HO densities give < >i/2radu of 
2.67 and 2.74 fm for ^^O and ^*0, respectively, which are 
rather close to those given by the IPM densities. Predic- 
tions given by the microscopic HFB densities are also 
in a good agreement with the data for even O isotopes 
excepting the ^''O case, where the HFB density gives ob- 
viously a too small < >^/^radius. Since the HFB cal- 
culation already takes into account the continuum effects 
[33], such a deficiency might be due to the static defor- 
mation of ^^O. A jump in the < >^/^value was found 
for ^■'O which could indicate to a neutron halo in this 
isotope. Behavior of < >^/^radii given by the best- 
fit HO densities agrees with the trend established by the 
OL of Glauber model [1] but, like the case of G isotopes, 
they are about 0.1 fm larger than those deduced from the 
OL of Glauber model. Thus, the OL of Glauber model 
seems to consistently overestimate ctr for the neutron- 
rich C and O isotopes under study in comparison with 
our approach. 

One clear reason for the difference between our results 
and those given by the OL of Glauber model analysis is 
that one has matched directly the calculated ctr with the 
measured ai in the Glauber model analysis [1] to deduce 
the nuclear radius. If we proceed the same way with 
the HO densities for the considered nuclei, the best-fit 
< >^/^radii decrease slightly but are still larger than 
those given by the OL of Glauber model. As already 
discussed in Sect. II, the zero- angle approximation for 
the NN scattering amplitude used in the Glauber model 
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might reduce significantly the strength of the exchange 
part of the imaginary OP given by Eq. (12) and could 
overestimate, therefore, the absorption in the dinuclear 
system. This effect should be much stronger if one uses a 
realistic finite-range representation of the NN scattering 
amplitude. Bertsch et al. have shown [30] that the zero- 
range approximation for the NN scattering amplitude 
leads to a reduction of the calculated ctr or an enhance- 
ment of the nuclear radius by a few percent (see Figs. 2 
and 3 in Ref. [30]). Thanks to such a cancellation of the 
exchange effects by the zero-range approximation for NN 
scattering amplitude, the simple OL of Glauber model 
was able to deliver reasonable estimates of nuclear radii 
for many stable and unstable isotopes [1]. It should be 
noted that the eikonal approximation for the scattering 
wave function used in the Glauber model was introduced 
in the past to avoid large numerical calculations. With 
the computing power available today, there is no problem 
to perform the OM and DWIA calculations for different 
nucleus-nucleus systems involving large numbers of par- 
tial waves, and the folding -|- DWIA method presented 
here can be recommended as a reliable microscopic ap- 
proach to predict the elastic scattering cross section and 
to deduce the nuclear radius from the measured ui. 

IV. CONCLUSION 

In this work we have explored the reliability of the op- 
tical model -h DWIA approach as a tool for extracting 
important information on nuclear sizes from interaction 
cross section measurements. We concentrate on the en- 
ergy region of 0.8 to 1 GeV/nucleon where interaction 
cross section data exist for various combinations of sta- 
ble as well as unstable projectiles on different targets. 
At these bombarding energies our knowledge of the em- 
pirical optical potential is scarce, especially for unstable 
systems, and we have used, therefore, the folding model 
to calculate the microscopic (complex) optical potential 
and inelastic form factors necessary for our analysis. 

We have chosen for the folding input the fully finite- 
range t-matrix interaction developed by Franey and Love 
[13]. The folded optical potentials and inelastic form fac- 



tors are used further as inputs for the standard optical 
model and DWIA calculations of total reaction cross sec- 
tions and interaction cross sections induced by stable and 
unstable He, Li, C, and O isotopes on ^^C target. By us- 
ing the well tested nuclear g.s. densities for the stable 
"^He, ^^C and ^^O isotopes, we found that the Franey 
and Love t-matrix gives extremely good account of the 
measured cti for these nuclei. 

We have further used the nuclear g.s. densities ob- 
tained in various structure models to calculate aj and 
made realistic estimate for the nuclear radii of (still 
poorly known) neutron-rich isotopes based on the com- 
parison between af^^'^ and a^^^. For the chains of C and 
O isotopes, our results agree reasonably well with the 
empirical trend established by the OL of Glauber model 
[1], but give consistently larger < >^/^radii for these 
nuclei. Such an effect could be due to the unsatisfactory 
treatment of the exchange part of the nucleus- nucleus OP 
in the Glauber model calculation. 

Although the nuclear radii deduced by our approach 
for some light halo nuclei might be a few percent smaller 
than realistic values because the dynamic few-body cor- 
relation was not considered explicitly in the present fold- 
ing -|- DWIA formalism, this fully microscopic approach 
was shown to be more accurate than the OL of Glauber 
model. Given realistic nuclear densities, it can give a re- 
liable (parameter-free) prediction of the nucleus- nucleus 
optical potential at energies around 1 GeV/nucleon. 
Therefore, it provides the necessary link to relate the 
calculated ai to the nuclear density and RMS radius. 
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